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Abstract 

We study the extension of our translationally invariant treatment of few-body nu- 
clear systems to heavier nuclei. At the same time we also introduce state-dependent 
correlation operators. Our techniques are tailored to those nuclei that can be dealt 
with in LS coupling, which includes all nuclei up to the shell closure at ^ = 40. 
We study mainly p-shell nuclei in this paper. A detailed comparison with other 
microscopic many-body approaches is made, using a variety of schematic nuclear 
interactions. It is shown that our methodology produces very good energies, and 
presumably also wave functions, for medium mass nuclei. 



1 Introduction 



During the last few years there has been an impressive improvement in the 
methodology used to perform fully microscopic calculations for the ground- 
state properties of nuclear systems. Most progress has been made for light 
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nuclei, particularly the systems of three and four nuclcons, where quite dif- 
ferent methods, including Monte Carlo methods [1]. Faddeev and Faddeev- 
Yakubovski methods [2], hyperspherical harmonic expansions [3], and transla- 
tionally invariant configuration interaction methods [4], have been pushed to 
a high level of rehabihty. 

Ideally one would like to extend as many of the methods as possible to heavier 
systems. Unfortunately, however, some of these methods are specially designed 
for very light systems, and cannot easily be extended to heavier systems. 
The Green's function Monte Carlo methods have been extended to systems 
of up to six nucleons, but it is unlikely that such methods will be able to 
tackle systems with more than a few additional nucleons in the foreseeable 
future. The hyperspherical harmonic method is also limited in the number of 
nucleons it can treat, even though the Pisa group has started an ambitious 
project for p-shell nuclei. Arias de Saavedra et al. [5] have recently investigated 
the application of the Fermi hypernetted chain (FHNC) method to heavier 
nuclei, with some apparent success. Another promising approach is one based 
on a stochastic variational principle [6], but this method is also limited to 
relatively light nuclei. We should also mention the variational Monte Carlo 
calculations of ^^O [7,8] which have considered reahstic interactions. These 
impressive calculations, however, have used a multiplicative cluster expansion 
to deal with the operatorial structure of the trial wave function. 

In this paper we shall discuss another approach, which was pioneered by us 
some years ago [4]. This translationally invariant configuration interaction 
(TlCl) method was inspired by the coupled cluster (CC) method. In its lowest- 
order implementation (CC2) the coupled cluster method consists of using wave 
functions obtained by acting on a reference state with the exponentiated one- 
and two-body correlation operators, expfS"! -|- S2\. The CC method is most 
naturally formulated in the occupation-number representation. However, as 
was to be expected, we found very slow convergence with respect to the cut- 
off on the single-particle basis in this representation. This led us to consider 
a linearised form of the CC2 approach, namely the TICI2 method, which also 
includes pair correlation effects but which can easily be converted to coordinate 
representation. It has the trial wave function 



In order to remove spurious effects due to the centre-of-mass motion the uncor- 
related wave function l^^o) must be a harmonic-oscillator shell-model function. 
In contrast to the standard CC2 method, the TICI2 method is variational, and 
leads to an upper bound to the ground-state energy. 

The simple form of eq. (1) for the wave function has provided very good es- 




(1) 
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timates of the ground-state energy of "^Hc for simple Wigner forces [4]. When 
we revert to occupation- number representation, we find that the Ansatz (1) 
corresponds to the most general shell-model excitation operators containing 
both one particle-one hole and two particle-two hole components. These com- 
ponents have to be appropriately related so as to ensure the translational 
invariance. 

One obvious improvement of this Ansatz is the introduction of three-body 
correlations, of the form 

I*) = (E/(^*j)+ E 9{n„Uk,rjk)] l^o). (2) 

\i<j i<j<k J 

Our calculations using this form of trial wave function, with a simple model 
interaction in "^He, produced almost the exact energy for the ground state [4]. 
As one can see, we have a systematic way of improving the wave function, 
by adding terms of physical relevance. One may expect that, in medium light 
nuclei, one must include only up to three- and maybe four-body correlations, 
in order to get a precise determination of the energy and the structure of the 
system. 

On the other hand, one must also deal with the complex operatorial structure 
of realistic nucleon-nucleon (NN) interactions, which in turn requires an ad- 
ditional operatorial structure for the correlation operator. At the same time, 
one also needs to find a way to extend the method to deal with nuclei with a 
larger number of particles. As is also true for some (but not all) of the other 
approaches cited above, this extension is trivial, but the underlying computa- 
tional complexity requires an appreciable effort. We should recall the obvious 
fact that the expectation value of a two-particle operator, like the NN po- 
tential, between states of the type of eq. (1) requires the evaluation of the 
expectation value of up to six-body operators. These calculations are large 
but not prohibitively so. 

The dual purpose of this work is to show how this kind of computation may 
be carried out in practice, and also to investigate the importance of the op- 
eratorial structure of the pair correlation. In this paper we are considering 
simple interactions of the V4 form (central, but spin- and isospin-dependent) . 
We also consider a two-particle correlation operator of the same structure, i.e., 
with the spin- and isospin-dependent form 

fij = fMj) + fa{rij){Oi ■ <?j) + fr{rij){T, ■Tj)+ faT{nj){ai ■ aj){fi ■ fj).{3) 

In sec. 2 we discuss the computational problems associated with this form 
of correlation operator and our way of tackling them. Our scheme goes be- 



3 



yond the typical shell-model methodology, which basically deals with two- 
body operators. In its current form our calculations are limited to selected 
nuclei, specifically to those systems which may be described or approximated 
by spin-isospin saturated wave functions. At the cost of losing the rotational 
invariance, one may basically deal with A — An nuclei, where n is an integer. 
For some systems (like, e.g., ^^C or ®Be) it is actually better to break the rota- 
tional symmetry in favour of a deformed structure. In these nuclei the energy 
of the deformed state is much lower than that built on a spherical state of 
the jj coupling scheme. A further restriction on the nuclei we can deal with 
originates from the fact that our method makes intricate use of L5'-coupling. 
This means that our method will not perform well beyond the shell closure at 
40Ca. 

Having discussed our methodology in some considerable detail, we then turn 
to a comparison of results in sec. 3. Our main goal is to gain some insight into 

both the power of our method and any associated shortcomings. We make a 
detailed comparison with other many-body methods and exact results where 
available. Finally we draw some conclusions and give an outlook in sec. 4. 



2 The evaluation of matrix elements 

In this section we shall show how to calculate matrix elements of a V4-type 
potential, when the correlations are of the form given in eq. (3). We present all 
reductions of the required 74-body matrix elements by identifying the minimal 
set of integrals that we are required to calculate. The reduction of spin and 
isospin degrees of freedom is also discussed, for the case of saturated systems. 

The discussion has been limited to the case of the expectation value of the 
potential, which is the most complex quantity we need to evaluate. A similar, 
but much less complicated, analysis can be carried out for the kinetic energy 
operator as well as for the computation of the norm of the pair-correlated 
state. 



2.1 Diagrammatic analysis 

Let l^o) represent an independent-particle harmonic-oscillator shell-model 
wave function for the nucleus under consideration. For the case of the spin- 
isospin saturated systems in which we are interested, this uncorrelated state 
may be represented by a single Slater determinant, which is hence fully char- 
acterised by the enumeration of the occupied single-particle states. Moreover, 
given that we are dealing with spin- and isospin-saturated systems, it is only 
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necessary to specify the harmonic oscillator states occupied, and we can ignore 
spin and isospin for the characterisation of the states. For example, the uncor- 
related state for ^^O requires the four single-particle states Os, Op+i, Opo and 
Op_i, where the subscript refers to the magnetic quantum number. We can 
transform to a Cartesian basis where the occupied single-particle states are 
represented by the three harmonic oscillator quantum numbers {nx,ny,nz)- 
The reference state for ^^O in this basis consists of Slater determinant formed 
from the four occupied states (0, 0, 0), (1, 0, 0), (0, 1, 0) and (0, 0, 1). 

Consider now the expectation value of the potential energy operator, given by 

{V) = {^o\J:fL{^J)J:V{kl) Y: fR{mn)\^o), (4) 

i<j k<l rn<n 

where the subscripts L and R identify (for later ease of discussion) the corre- 
lation operators at the left hand side and the right hand side of the potential 
operator, respectively. The indices in parentheses label the particles involved in 
the corresponding correlation or potential operator. When we expand the sums 

in eq. (4) there appears a total of ( ^ ) terms. One way to classify the individ- 
ual terms is by the number N of particles involved. This number ranges from 
2 (e.g., a term A(12)y(12)/^(12)) up to six (e.g., a term /^(12)y(34)/«(56)). 
Moreover, in order to reduce the A-particle matrix element to an A^-particle 
matrix element it is convenient to decompose the operator of eq. (4) into a 
sum of symmetric A/'-body operators. It is very simple to make a diagrammatic 
analysis of these operators by assigning to all /'s and V^s a hue connecting 
the points specified in their arguments [4] . The topologically distinct diagrams 
constitute a natural basis for the expansion of the matrix elements. 

With some patience one then obtains the decomposition 

16 

E/4U)E^(^0 E fnirnn) = E E ^^n....... (5) 

i<j k<l m<n d=l ii<i2<— <2jv^ 

where the symmetric operators are defined in table 1. To clarify the mean- 
ing of this table, let us evaluate, for example, the three body part of JlV/r- 
According to the table it can be expanded in the diagrams 2 till 5 as 

(/LV^/ii)iii2i3 =E[/i(^i^2)V^(«i«2)/ii(«i«3) + fL{.Hi2)V{iiii)fR{iiiz) + 
p 

fL{iii2)V{iii3)fR{iii2) + fL{iii2)V{iii3)fR{i2i3)] 

where the sum over permutations P includes the six permutations of the triplet 
(^i^2^3)- In other cases there may be restrictions on the permutations, as men- 
tioned in the table. The number of permutations of a given diagram under 
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these restriction is listed in the last column of the table. This allows the enu- 
meration of the diagrams to be checked, since they must satisfy the equality 





where d refers to a given diagram, and is the number of particles involved 
in that diagram. 

After this classification, the expectation value of a given diagram between 
Slater determinantal wave functions is reduced to 



ii<...<iN^ 

= E (V^n(l)---^^.,TO|0LivjEeP^Pn(l)---^P«.,TO).(6) 

^^d- ji...ijv^=l P 

This equation requires further explanation. The single-particle states are rep- 
resented by the wave functions ipi{j), which also includes a spinor and an 
isospinor. Here the subscript i refers to the single particle state, and the index 
in parenthesis specifies the particle involved. The right-hand side of the equa- 
tion also includes a sum over permutations of the single-particle state labels, 
without any restriction. Note that the sum over single particle states extends 
over all occupied orbitals. 

A very important reduction in the computation is related to a further sym- 
metry of eq. (6), and refers to the sum of permutations (not explicitly written 
in this equation) which enter into the definition of O'^. It turns out that each 
of these permutations gives the same contribution to the matrix element, so 
that it is not necessary to evaluate all permutations separately. It is sufficient 
to consider one of the diagrams contributing to O'^, and multiply the resulting 
expectation value by the number of permutations of the diagram given in the 
last column of table 1. 



2.2 The spin and isospin degrees of freedom 



To evaluate a given matrix element, eq. (6), one must sum over all occupied 
single-particle states. In the case of spin-isospin saturated systems, the index 
labelling the single-particle states may be split into three pieces, 

^=(^^^^^-) 
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Table 1 

The topologically distinct diagrams for potential matrix elements. The column la- 
belled d gives a unique label for each diagram, specifies the number of parti- 
cles involved. The column labelled fi gives the particle labels (ij) of the left-hand 
side correlation function, V those (kl) of the potential and /r those (mn) of the 
right-hand side correlation. The diagrams must be symmetrised with respect to the 
particle labels ii - . involved, but with some restrictions to avoid overcounting. 
These restrictions are given by the order in which the labels appear on the functions 
fL,R and V. The number of permutations Pd under these restrictions is given in the 
last column. 



d 


Nd 


h 


V 


fR 


Restrictions 


Pd 






ij 


kl 


mn 


on Permutations 




1 


2 


i\i2 


kk 


kh 


i < j 


1 


2 


3 


i^ in 


k io 


i^ i'^ 
i "o 


None 


6 


3 






kh 


kh 


None 


6 


4 




hh 


kh 


kk 


None 


6 


5 




i\i2 


kis 


kh 


None 


6 


« 

D 


A 


i\i2 


kk 


hu 




a 
o 


7 




hh 


hu 


hu 


i < j, m < n 


6 


8 




ki2 


hu 


kk 


i < j, k < I 


6 


9 




ki2 


kh 


hu 


None 


24 


10 




ki2 


kh 


ku 


None 


24 


11 




ki2 


hu 


kh 


None 


24 


12 




ki2 


kh 


ku 


None 


24 


13 


5 


ki2 


hu 


kh 


k < I 


60 


14 




ki2 


kh 


iih 


m <n 


60 


15 




ki2 


hu 


hh 


i < j 


60 


16 


6 


ki2 


hu 


hh 


i < j, k < I, m < n 


90 



where s, a and r refer to space, spin and isospin, respectively. The sum over all 
occupied states is equivalent to a sum over harmonic-oscillator states, together 
with a replacement of all products of spin and isospin operators by their 
respective traces. 

The computation of the trace requires a further reduction, since it interferes 
with the calculation of permutations. Each of the permutations appearing in 
the A'"-particle state of eq. (6) can be split into three independent permuta- 
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tions, as can be seen from the pair-wise exchange 

The pairwise permutations referring to spin and isospin indices can be imple- 
mented by having famihar spin and isospin exchange operators, P[j and P[j 
act on the wave function. Thus the spin and isospin part of the permutation 
P, P"^ and P'^, can be written as products of exchange operators, which can 
be included inside the trace over spin and isospin. 

We thus obtain the matrix element 
(*o|/L(^J)^(A;/)/i^(mn)|^o) 

A/4 

pqr j^...t=^,=l P 

(0,5(1) . . .0.^(iV)|/£(r,,-)^^(r«)/fl(r.„n)|0pq(l) ■ ■ .<l>p^l,{N)). 

Here the 0's are the spatial (harmonic oscillator) parts of the single-particle 
wave functions ■0 in eq. (6). Since all spin and isospin sub-states are filled, each 
harmonic oscillator state can be occupied at most four times. Both the left 
and right correlation operators have been expanded in a basis of operators, 

fi^j) = T.f'in,)Q'{tj), (8) 

k 

with a similar expansion for the potential operator. A quite convenient basis 
for a V4 form is the set of spin and isospin exchange operators, 

e\tj) = 1, Q\zj) = P^^, e'izj) = P[^, Q\zj) = P^P^. (9) 

Since the spin and isospin permutation operators can also be written as prod- 
uct of exchange operators, we now have to evaluate the trace of a long string 
of exchange operators. 

2.3 The computation of traces 

For the case of a V4 set of operators, both for correlations and interaction, 
the calculation of the traces is rather simple. For an A^-particle subsystem, 
a given spin (or isospin) state may be represented by a column vector of 2^ 
components, where a given component defines all orientations of the N spins. 
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A simple representation is obtained by assigning a binary sequence to a given 
set [(Ti, (72, . . . , (7iv] of spins: 

= (71 + 2(72 + 4(73 ... + 2(^-^)(7Ar (10) 

where we use ai — for spin up, and (7, = 1 for spin down. 

With this notation, a spin exchange operator is represented by a 2^ x 2^ ma- 
trix with only one non-zero element per row and column. This facilitates the 
storage of matrices, the multiplication of exchange operators and the compu- 
tation of the trace. 

We have to evaluate a large number of traces, since we have four compo- 
nents for the left correlation, the potential and the right correlation, as well 
as the 720 permutations of a six-particle subsystem, leading to a total of 
720 X 4^ = 46080. With the representation introduced above it is however 
rather straightforward - and very efficient - to write a computer program 
that evaluates all these traces. 

2.4 Spatial integrals 

In order to make the calculation manageable we expand the correlation func- 
tions in a set of Gaussians, 

/"(r)=Ec?exp(-Ar^), (11) 

i 

a technique that has been proven to give an excellent representation of the 
correlation function / if we include negative as well as positive values of (3 [4]. 
This is well matched to the use of harmonic oscillator single-particle states, 
since the required space integrals can be shown to be of the form of the product 
of an exponential of a positive-definite quadratic form with a polynomial. 

If the potential is a combination of Gaussians, hke the S3 interaction of Afnan 
and Tang [9], it may also be absorbed into the quadratic form in the expo- 
nential. Then, the required integrals can be computed by using a recurrence 
relation [10]. 

For other algebraic forms, we can follow two routes. One is to fit a set of Gaus- 
sians to the potential, and apply the technique discussed above. We choose 
to follow another route [10], where we perform all integrals apart from those 
involving the coordinates in the potential by the technique sketched above. 
The remaining one-dimensional integral may then be computed by means of 
a suitable numerical method. 
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2. 5 Solution 



Since the different components in the decomposition of eq. (11) of the corre- 
lation function are not orthogonal, the expectation value of the energy can be 
given in the form 



E{{4}) = '"^"^ (12) 



(jp),(i9) 



where cf is the coefficient multiplying the ith Gaussian in the correlation 
function attached to operator p. If we now optimise with respect to cf * we end 
up with a generalised eigenvalue problem 

E^IT^ = ^E-A/^M- (13) 

{ji) in) 

The lowest eigenvalue A gives the approximation for the ground-state energy. 
2.6 The helium case 



The four-nucleon system deserves a special treatment, since the decomposition 
of the two-particle correlation introduced in eq. (8) is over-complete for this 
special nucleus. This fact is related to the special form of the uncorrelated 
state, which is fully space-symmetric, so that the action of the operators tied 
to the correlation function is strongly simplified. In particular, the action of 
the spin-isospin exchange operator PijP[j is equivalent to —1, and the action 
of P^j is the same as the action of —Plj- Clearly, one should limit the general 
form of the correlation so as not to create duplicate basis states, since this 
will give rise to problems when solving the generalised eigenvalue problem of 
eq. (13). Specifically, the state-dependent correlation of the V4 type for ^He 
should include only the central scalar (Wigner) and spin-exchange (Bartlett) 
pieces. Equivalently, one may include two separate correlation functions, one 
for the singlet channels and one for the triplet channels. 



3 Results 



Using our technology we have set out to calculate the binding energy of *^He, 
^Be, ^^C and ^^O. The results are presented in table 2. In order to allow for 
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Table 2 

Energies of p-shell nuclei for the five interactions considered, Bl, S3, MS3, MT 
I/III and MT V. The column labelled SI corresponds to a purely scalar (state- 
independent) pair correlation, and the column labelled SD corresponds to a (state- 
dependent) pair correlation of V4 type. 



Bl interaction 


Nucleus 


a (fm 


-1) Esi (MeV) 


^SD (MeV) 


^Hc 


0.729 


-37.86 


-37.86 


SBe 


0.595 


-49.18 


-61.30 


12c 


0.595 


-84.91 


-103.93 


160 


0.602 


-145.94 


-167.30 


S3 interaction 


^He 


0.717 


-25.41 


-28.19 


SRe 


0.633 


-34.15 


-44.26 


12c 


0.671 


-71.14 


-87.78 


160 


0.707 


-141.64 


-164.88 


MS3 interaction 


^He 


0.713 


-25.41 


-27.99 


SBc 


0.582 


-26.26 


-37.30 


12c 


0.588 


-46.22 


-62.99 


160 


0.596 


-85.56 


-105.64 


MT I/III interaction 


^He 


0.741 


-29.45 


-30.81 


8Be 


0.660 


-46.67 


-52.67 


12c 


0.698 


-99.05 


-109.04 


160 


0.744 


-194.10 


-207.52 


MT V interaction 


^He 


0.741 


-29.45 


-29.45 


SBe 


0.869 


-129.25 


-130.23 


12c 


0.997 


-425.99 


-429.44 


160 


1.078 


-966.65 


-973.67 



the best possible comparison with existing results in the literature we have 
investigated the Brink-Boeker Bl potential [11], the Afnan-Tang S3 potential 
[9], the modified S3 (MS3) potential [12], and the Malfiiet-Tjon MT I/III and 
MT V potentials [13]. We have optimised the value of a, the inverse length- 
scale of the harmonic oscillator, for the case of the state- independent (SI), 
purely central (scalar) correlations, i.e., where we restrict the expansion of 
eq. (8) to the k = 1 term only. We then performed a calculation with the 
full state-dependent (SD) V4 type correlation operator, keeping the value of 
a fixed. The difference between the results from the use of a central state- 
independent (SI) and the use of a state-dependent (SD) correlation operator 
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is thus probably slightly underestimated. Our calculations for ^Bc and ^^C 
were performed using a non-spherical reference state (in Cartesian notation) 



v|>o) = (0,0,0)^(0,0,1)^ 



(14) 



for ^Be and a similar state, 



(0,0, 0)^(1,0, 0)^(0, 1,0)^ (15) 



for ^^C. We recall that wc can only use states of this form since we work 
so as to have spin and isospin saturation, within an LS coupling scheme. We 
have broken the rotation symmetry in these cases, and we no longer get a strict 
estimate for the ground-state energy, but for the slightly higher intrinsic energy 
- a weighted average of the energies of the J = 0, 2, 4 states comprising the 
ground-state band. Nevertheless, such an estimate is preferable to one using 
a spherical wave function (e.g., (0^1/2)^(0^3/2)^^ for the ground state of ^^C). 
This can already be seen from a naive calculation of the ^^C nucleus using 
harmonic oscillator states and the Bl potential, without correlations. The 
energies in the spherical state (0si/2)''(0]93/2)®^ and in the non-spherical state 
I (0, 0, 0)^(1, 0, 0)^(0, 1, 0)^) are -43.11 MeV and -59.04 MeV, respectively, while 
a shell-model (or configuration interaction) diagonalisation of the Hamiltonian 
in the Op-shell gives -62.07 MeV. Thus an uncorrelatcd estimate based on the 
deformed state appears to be almost optimal. Note that the energy may be 
lowered still further by using a deformed single particle basis, i.e., by using 
different harmonic oscillator parameters along the different axes. 

The effects of improving the correlations are very interesting. Note that for the 
Bl and MT V interactions the introduction of SD-type correlations does not 
make a difference for ^He. This is due to the fact that these interactions contain 
only purely central scalar (Wigner) and space-exchange (Majorana) terms, and 
hence do not couple to the spin-dependent piece of the correlation operator 
in "^He. Even though there is no spin-dependent part to these potentials, spin- 
dependent correlations do play a role for all heavier nuclei, so long as the 
uncorrelated state is not fully space-symmetric. In the Bl case the difference 
is over 20 MeV for ^^O. The effect is smaller for the MT V potential, but this 
interaction is very pathological - it does not saturate nuclear matter. For the 
S3, MS3, and MT 1/111 interactions we find that the effect of the additional 
correlations is in the 10 to 20 MeV range for ^^O, an important difference. 

In figure 1 we show the correlation functions for the four p-shell nuclei. Apart 
from the case of ®Be the central scalar correlation function changes little from 
the SI to SD calculation. Indeed, one could almost keep this one fixed and 
only vary the additional three. This may be due to some kind of perturbation 
argument: in general the spin- and/or isospin-dependent correlation functions 
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1.0 



1.0 



012345 012345 

r(fm) r(fm) 

Fig. 1. The calculated correlation functions in the case of the MS3 potential. Dashed 
lines correspond to purely central scalar SI correlations and solid lines to SD cor- 
relations of V4 type. In the latter case the labels 1, a, r and err correspond to the 
terms A; = 1,2, 3 and 4 in the parametrisation of eqs. (8) and (9). As discussed in 
the text, for helium only two such functions should be used. 



arc a lot smaller than the central one. The spin and isospin correlation func- 
tions are almost equal in magnitude but opposite in sign. There is no obvious 
argument why this should be so, but it may be due to some approximate 
persistence of the corresponding exact relation in the case of ^He. 

In order to be able to judge the power of our technology we compare sev- 
eral different methodologies for the purely scalar MT V potential in table 3. 
For ^He our results compare favourably with the Green's function and diffu- 
sion Monte Carlo results, both of which have now been superseded by various 
highly accurate calculations, including one based on the stochastic variational 
principle (see ref. [6] and references in that paper). More surprising is the com- 
parison for ^^O. For the MT V potential, with its extremely strong binding, 
^^O is a high-density system. Our calculations with only two-body correla- 
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Table 3 

A comparison of the binding energies (in MeV) of ^He and ^^O using various tech- 
niques with the Malfiiet-Tjon MT V potential. TICI2 is the method presented in 
this paper, VMC are various variational Monte Carlo calculations, GFMC and DMC 
are fixed-node Green's function and diffusion Monte Carlo methods, FHNC is the 
Fermi hypernetted chain method based on a Jastrow wave function, IDEA is the so- 
called integro-differential equation approach, and SVM is the stochastic variational 
method. 





^He 


16q 


TICI2(SI) 


29.45 


966.65 


TICI2(SD) 


29.45 


973.67 


VMC 




1024 ± 5 [14] 






1103 ± 1 [15] 






1138.5 ± 0.2 [16] 


GFMC [14] 


31.3 ±0.2 


1194 ± 20 


DMC 


31.32 ± 0.02 [17] 


1189 ± 1 [16] 


FHNC/0 [18] 




987 / 1152 


FHNC//0 [19] 




1059 / 1055 


IDEA [20] 


30.7-31.2 


1021-1027 


SVM [6] 


31.360 





tions are able to reproduce over 80% of the binding energy (974 MeV for 
TICI2(SD) as compared to 1189 MeV for the very accurate DMC result). 
Inclusion of state-independent three-body correlations in our calculation - a 
relatively straightforward extension - should probably be able to give a highly 
accurate answer. 

The TICK method also appears to be competitive with several of the methods 
based on Jastrow-correlated non-interacting wave functions of a single Slater 
determinant form, where the Jastrow correlation operator is a product over 
all pairs of particles of a scalar (state-independent) pair function of the in- 
terparticle distance only. The many-body energy expectation value for such 
trial Jastrow wave functions cannot be calculated analytically, but may be cal- 
culated within controlled error bands using variational Monte Carlo (VMC) 
techniques. Each of the three VMC results cited in table 3 is of this form, with 
the differences between them resulting from different choices of single-particle 
wave functions and from different forms of the Jastrow correlation function 
over which the variational search is made. These are discussed in more de- 
tail below. As an alternative to the stochastically exact VMC calculations of 
the energy expectation value, the various Fermi hypernetted chain (FHNC) 
approximations [5,18,19] represent massively resummed cluster expansions of 
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the same quantity, in which some cluster terms (of chain or ladder form) arc 
summed to all orders but various so-called elementary diagrams are neglected 
or approximated. The resulting FHNC approximants for the ground-state en- 
ergy are no longer strict variational bounds, insofar as the higher-order terms 
have been neglected. 

Firstly, the VMC resuh of 1024 ± 5 MeV for the binding of ^^O with the 
MT V potential from ref. [14] uses single-particle wave functions generated by 
a Woods-Saxon potential and a scalar Jastrow correlation function obtained 
from a constrained Euler- Lagrange optimisation of the second order cluster 
expansion approximation to the energy functional. We note that the corre- 
sponding result from the lowest order FHNC calculation, in the formulation 
of Fantoni and Rosati, the so-called FHNC/0 method [18] in which the el- 
ementary diagrams are neglected completely, using the same single-particle 
wave functions and Jastrow correlation function, is 987 MeV. Secondly, the 
VMC result of 1103 ± 1 MeV from ref. [15] for the ^^O binding energy em- 
ploys harmonic oscillator single-particle wave functions and a Jastrow correla- 
tion function of Gaussian form, /(r) = 1 -|- aexp(— with the oscillator 
length parameter and the constants a and b optimised. Finally, the similar 
VMC result of 1138.5 ± 0.2 from ref. [16] uses optimised harmonic oscillator 
single-particle wave functions and a Jastrow correlation function of the form 
/(r) = exp(-ae"''^). 

We note that a so-called FHNC/ /O calculation for ^^O using the MT V poten- 
tial has been performed by Krotscheck [19]. By using the best possible choices 
for both the single particle wave functions and the Jastrow correlation func- 
tion (i.e., by solving the corresponding Euler-Lagrange equation which arise 
by making the FHNC/ /O energy functional stationary) he obtains a value of 
1059 MeV for the binding energy. Furthermore, by optimising only the Jastrow 
correlation function via his Euler-Lagrange approach, and using single-particle 
wave functions obtained from a Woods-Saxon potential which reproduces the 
r.m.s. radius of the GFMC calculation of ref. [14], he obtains the only very 
slightly worse value of 1055 MeV. He estimates, using the FHNC / C approxi- 
mation that the contribution from the neglected exchange diagrams is 14 MeV. 
We also note that by employing precisely the same nearly optimal Woods- 
Saxon single-particle wave functions, but with a Jastrow correlation function 
obtained from a constrained minimisation of the two-body cluster expansion 
approximation to the energy functional, a corresponding FHNC/0 calculation 
[18] gives a ^^O binding energy of 1152 MeV. This difference of nearly 100 
MeV is presumably a measure of the uncertainty in the FHNC/0 results, in 
the light of the other results cited. It seems rather unlikely that the corre- 
sponding VMC calculation for this same wave function will give such a large 
value for the binding energy. It would thus appear that the latter FHNC/0 
result does not provide a true variational estimate. We defer to sec. 4 a further 
discussion of the FHNC results. 
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Table 4 

A comparison of the binding energies (in MeV) of ^^C and ^^O using various tech- 
niques with the Brink-Boeker Bl potential. For a description of the methods see 
the caption of the previous table. FAHT is a cluster expansion. 





12 c 


16 


TICI2(SI) 


84.91 


145.94 


TICI2(SD) 


103.93 


167.30 


VMC [21] 


82.9 ±0.2 


150.9 ±0.3 


FHNC/0 




168.2'^ / 167.2'' [18] 




68.9'=[5] 


161.7 [5] 


FHNC-1 




150.4'^ / 152.4'' [18] 




63.7^ [5] 


151. 4^= [5] 


FAHT-4 [21] 


82.3 


148.6 


IDEA [20] 


79.8-80.1 


164.4-165.2 



" state-independent Gaussian, LS coupling 



statc-indcpcndcnt 2-body Eulcr-Lagrange, LS coupling 
state-independent 2-body Euler-Lagrange, jj coupling 

The reason for investigating the MT V potential is that quite a few results are 
available, including very good GFMC and DMC calculations. However, we also 
investigate the Bl and MS3 potentials which have more sensible saturation 
properties, but for which no "exact" results are available to act as a benchmark 
for the calculations. Bl is a Wigner-plus-Majorana potential, whereas MS3 is 
a full V4-type potential. It is interesting to look at both potentials in turn, 
since they highlight different aspects of the use of V4-type correlations. 

For the Brink-Boeker Bl potential we present results in table 4 for the nuclei 
^^C and ^^O, the only nuclei for which several other results are available. For 
^^C we must distinguish between calculations like our own that studied the 
intrinsic state for the ground-state band and those that start from a spherical 
state. Unfortunately, there are no GFMC or DMC results available for either 
nucleus with the Bl potential, and we can only compare with the results from 
other many-body techniques. All of the other results in table 4, except those 
from IDEA, are based on a Jastrow form of the wave function, again using 
a scalar (state-independent) correlation operator. The results labelled FAHT- 
n are from an nth order cluster expansion of the Jastrow energy functional 
in the so-called factorised Aviles-Hartogh-Tolhoek decomposition [22], which 
is multiplicative van Kampen-like version of the original additive Ursell-like 
AHT version. 
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Our SI calculation for ^^C results in slightly better values for the binding 
energy than the VMC and FAHT calculations. Note, however, that Jastrow 
calculations with anisotropic HO wave functions result in an increase of the 
binding energy of around 10 MeV (92.3 ± 0.3 MeV with the VMC method 
of [21]). However, when we include state-dependent (SD) correlations in our 
TICK calculation we find an increase in binding of nearly 20 MeV. Bearing 
in mind that our results are fully variational it is clear that our TICI2(SD) 
calculation is better than any of the other results shown for ^^C with the Bl 
potential. We note also that the results based on a spherical assumption give 
appreciably less binding for ^^C. This is particularly so for the FHNC results 
[5] based on a jj coupling scheme. We show results for both the lowest-order 
FHNC/0 approximation of the Fantoni-Rosati scheme, in which all elementary 
diagrams are neglected, and for their FHNC-1 approximation in which the 
first-order elementary exchange diagram has been included. Such diagrams 
are clearly likely to be of more importance for potentials of the Bl type with 
an appreciable Majorana component than for purely scalar potentials of the 
MT V type. The effect of their inclusion is seen to reduce the binding. The 
discrepancy between the VMC resuh [21] of 82.9 ± 0.2 MeV and the FHNC-1 
result [5] in the jj basis of 63.7 MeV is surprisingly large. We note, however, 
that the state-independent Jastrow calculations in ref. [5] were obtained in a 
calculation which included the Coulomb force, by minimisation of the energy 
at the second order of the cluster expansion. The results shown in table 4 have 
been obtained by subtracting out the quoted values of the Coulomb energy 
to the binding energies given in ref. [5]. Clearly, it is conceivable that had the 
minimisation been done without the inclusion of the Coulomb force, slightly 
improved results might have been obtained. We also note that the calculations 
in ref. [5] were performed with Woods-Saxon single-particle wave functions, 
thereby making it impossible to remove centre-of-mass effects exactly. 

In the case of ^^O with the Bl potential we find similar results to the ^^C case. 
Again our TICI2(SI) calculation gives a good result, and the improvement 
obtained by allowing for state-dependent correlations leads to what is possibly 
the best (and certainly the variationally most reliable) result now available 
for this case. The difference in binding between our SI and SD calculations 
of over 20 MeV is clearly very appreciable. We note that the FHNC/0 and 
FHNC-1 results for ^^O have been done in both LS [18] and jj [5] couphng 
schemes. For the former case we cite results using Jastrow correlation functions 
both parametrised as a Gaussian form, /(r) = 1 -|- aexp(— r^/6^), and from a 
constrained Euler-Lagrange minimisation of the second-order cluster energy. 
Woods-Saxon single-particle wave functions were used for all the FHNC results 
cited. 

We turn lastly to the case of the MS3 potential, and we again compare our 
results with others for this potential for the same two nuclei, ^^C and ^^O. 
We note that results for these nuclei for the MS3 potential are also available 
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Table 5 

A comparison of the binding energies (in MeV) of ^^C and ^^O using various tech- 
niques with the MS3 potential. For a description of methods see the caption of 
tables 3 and 4. BHF is the Brueckner-Hartree-Fock method. 



12c 

TICI2(SI) 46.22 85.56 

TICI2(SD) 62.99 105.64 

FHNC/0 113.5" [18] 

36.5754.2^= [5] 108.17145.5^= [5] 

FHNC-1 105.3" [18] 

34.7^ [5] 103.0'' [5] 

FAHT-3 (deformed) [12] 51.4 107.7 

BHF [12] 118.6 

IDEA [20] 44.2-44.4 102.9-103.2 



" state-independent 2-body Euler-Lagrange, LS coupling 

state-independent 2-body Euler-Lagrange, jj coupling 
isospin-dependent 2-body Euler-Lagrange, jj coupling 

within the FHNC framework using the jj coupling scheme and allowing, for 
the first time, some state-dependence in the Jastrow correlation operator. More 
explicitly, the authors of ref. [5] incorporate isospin-dependence by allowing 
different pairwise Jastrow functions for the proton-proton, neutron-neutron 
and proton-neutron pairs. Otherwise the FHNC results shown in table 5 are 
calculated as explained above for the Bl potential. We note again that the 
calculations of ref. [5] have included the Coulomb force. The results in table 5 
have been obtained by subtracting out the cited Coulomb energy in each case. 
Once again, for ^^C our results are better than all the others. The effect of 
the isospin-dependence in the FHNC calculation is comparable in magnitude 
to the effect of the state-dependence in our calculations. However, the FHNC 
results [5] are less bound for ^^C than our TICI2 results, presumably at least 
in part because of their use of a spherical reference state. 

For the case of ^^O nucleus with the MS3 potential, the effect of the in- 
corporation of the state-dependent correlations in our TICI2 calculations is 
again to increase the binding energy by about 20 MeV. Naively, one might 
have expected a somewhat larger increase than for the Bl potential, since 
the MS3 potential exploits the full V4 state dependence. Our own TICI2(SD) 
results are in this case close to most of the other results cited, with the excep- 
tion of the isospin-dependent FHNC/0 result, where the gain in binding over 
the corresponding state-independent case is nearly 40 MeV. In view of the 
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non- variational nature of both the FHNC and Brueckner-Hartree-Fock (BHF) 
results, it is difficult to draw more detailed conclusions. 



4 Discussion and outlook 

Wc have demonstrated that the TICK methodology provides a very reasonable 
starting point for the calculation of the binding energies of light-to-medium 
mass nuclei. The important state-dependent pairwise correlation effects can 
be efficiently incorporated without destroying the variational bound. By con- 
trast, most competing many-body methods suffer from various uncontrolled 
approximations, which often include the loss of a strict upper bound on the 
ground-state energy, even when the calculation starts from a variational frame- 
work. Such is often the case for FHNC and other more ambitious calculations 
performed within the correlated basis function (CBF) approach. Since such 
techniques provide perhaps the main competitor to our own approach outlined 
here, it is worth describing in some detail the source of the uncertainties that 
can typically arise in FHNC-type calculations. 

In a nutshell, the standard FHNC formulation of Fantoni and Rosati sums 
all of the Jastrow cluster-expansion terms that can readily be summed by 
one-dimensional integral equations (see ref. [5] and references cited therein for 
details). One drawback of this procedure is that certain cancellations between 
elementary exchange diagrams are thereby ignored. The resulting FHNC/n 
scheme then classifies the only omitted diagrams, the so-called elementary di- 
agrams, according to their number of vertices or points. Thus, the FHNC/n 
approximation includes elementary diagrams with up to n vertices. The re- 
sulting prescription is both clear-cut and leads to relatively simple integral 
equations which are vahd for small distances. Nevertheless, the long-range 
parts are more problematic. In particular, for infinite homogeneous systems 
one neither expects to get nor obtains the correct behaviour of the liquid 
structure function S{k) as k 0. Moreover this prescription does not allow, 
in principle, for correct optimisation. 

In contrast to the Fantoni- Rosati FHNC/n scheme, Krotscheck has proposed 
a different classification which maintains the Pauli principle at each step. One 
possibility among several to accomplish this is to classify the exchange dia- 
grams with respect to the number of internal lines rather than vertices. The re- 
sulting lowest-order approximant in this scheme is termed FHNC/ /O. In some 
well-defined sense it sums the union of the equivalents of the ring diagrams 
with a bare particle propagator and bosonic ladder terms. Effectively, for the 
case of infinite homogeneous systems, it sums correctly and self-consistently 
the most important diagrams in coordinate space when short-range correla- 
tions are important, and in momentum space when long-range correlations are 
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important. In turn this leads to an inconsistency between the pair distribution 
function, g{r), and the hquid structure function, S{k). 

Prom the above brief discussion it should be clear that the main difficulty and 
uncertainty in FHNC techniques is in the treatment of elementary exchange 
diagrams. We have already seen in sec. 3 the differences between the FHNC/0 
and FHNC-1 results, where the latter incorporates only the first-order elemen- 
tary exchange diagram from the FHNC/4 set. We also note in this context 
that the FHNC/C scheme of Krotscheck provides another means of incorpo- 
rating the omitted diagrams. Basically this scheme exploits the fact that one 
knows the values of these diagrams in the case of an infinite homogeneous 
system at momentum transfers k = and k = 2k-F, where k-p is the Fermi 
momentum. The FHNC/C scheme uses these exact properties to estimate 
the value of the omitted terms by a simple polynomial interpolation between 
these limits. Thus FHNC/C is a superset of FHNC/0, but with an approxi- 
mate treatment of what the latter omits. Such a correction of the FHNC/0 
equations is necessary for the optimisation. We have already commented in 
sec. 3 on the discrepancy between the FHNC/0 and FHNC//0 results in the 
case of ^^O with the MT V potential. It has been shown [19] that inclusion of 
the C class of diagrams has only a marginal influence. 

Apart from the above difficulties and uncertainties over the treatment of ele- 
mentary exchange diagrams, the incorporation of state-dependent correlations 
into FHNC treatments is also very fraught with difficulty. This arises essen- 
tially due to the non-commutativity of the correlation operators in a Jastrow- 
type product wave function when the scalar Jastrow function is replaced by 
an operatorial sum. The isospin-dependent correlations introduced in ref. [5] 
very recently, and reported in sec. 3, have led to workable FHNC calculations 
only because in this case the corresponding operators commute among them- 
selves. The introduction of more complicated operatorial forms (e.g., of the 
sort used by us) will not provide such a straightforward extension. Our own 
TICI2 method does not suffer from these drawbacks. 

We now turn briefly to the alternative IDEA approach, which is formulated 
purely in terms of the hyper-radius and the relative coordinate of the two par- 
ticles needed to describe the pair correlations incorporated in the method. In 
this method, as currently practised, several uncontrolled and often unjustified 
approximations are made. As a consequence the results are non-variational, 
and usually one has a range of values depending on which of many approxima- 
tions has been made, none of which is fundamentally preferred to another. A 
recent investigation [23] of a simple model suggests rather forcefully that even 
if no approximations are made, the inclusion of the hyper-radial excitations 
in a many-body wave function that lies at the heart of the IDEA method 
does not provide a very cost-effective improvement. This study suggest that it 
would be much more efficient to include pair (and triplet, etc.) correlations to 
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the best of one's ability, as done in the present work. Prom this viewpoint our 
own method of incorporating pair correlations is far more systematic than in 
the IDEA methodology. 

In conclusion, our method works well for interactions and correlations of V4 
form. Of course one needs a much more complicated interaction to describe 
phase shifts well: the latest Argonne V18 interaction has 18 components. Even 
though it may not be immediately obvious what kind of correlations are im- 
portant when one wishes to deal with such a complicated force, it is clear that 
the first step is to look at interactions and correlations of the V6 form, where 
one includes the tensor force and tensor correlations. Such an investigation is 
now under way, and it appears that the methodology sketched in Sec. 2 can 
be extended to deal with such a more realistic force. Actually when we have 
V6 type correlations, we may already pick up the major part of the spin-orbit 
force as well; at this level there is no easy match between the operatorial 
structure of the correlations and that of the potential. 

In the end we shall also have to include three-body correlations. Since the 
effect of such correlations is not expected to be very large, one might hope 
that inclusion of state-independent three-particle correlations is enough. This 
would still lead to very complicated calculations, but we believe that we can 
certainly deal with such a scheme. Indeed, we have already reported such 
results for '^He using purely central (scalar) potentials [4]. Results for heavier 
nuclei with more realistic potentials will be reported in the future. 
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